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Abstract: Using a recently developed effective field theory for the interactions of nucleons 
; j ■ at non-relativistic energies, we calculate non-perturbatively Coulomb corrections to proton-proton 

scattering. Including the dimension-eight derivative interaction in the PDS regularization scheme, 
we recover a modified form of the Blatt- Jackson relation between the scattering lengths. The 
effective range receives no corrections from the Coulomb interactions to this order. Also the case 
of scattering in channels where the Coulomb force is attractive, is considered. This is of importance 
for hadronic atoms. 

> 

m 
(N 

in 

m ; 1 Introduction 

O 

G\ 

Effective field theories are constructed to give a complete description of interacting par- 
| tides in terms of only the quantum fields which can be excited below a characteristic 

energy scale A. Degrees of freedom of higher energies are represented by an expansion 
of the Lagrangian in terms of local operators of increasing dimensions. For the nucleon 
system at energies below the pion mass m n the effective theory will thus involve only the 
nucleon field and derivatives thereof. It was first constructed by Weinberg]]]] and has been 
investigated by many others [g]. 

rS 

■ In order for the effective theory to be useful it must come with well-defined counting 

rules so that one knows which operators to include and which to disregard in a calculation. 
Only then can one obtain a systematic and reliable expansion in p/A for any scattering 
amplitude characterized by the momentum p. In the original formulation of the effec- 
tive theory for nucleons one were faced with problems in this connection because of the 
unnaturally large S'-wave scattering lengths. This brings a smaller energy scale into the 
system and the standard regularization schemes needed to handle the divergences from 
loop integrations had difficulties reproducing consistent counting rules |J. A solution of 
these problems was subsequently given by Kaplan, Savage and Wise with the introduction 
of a new regularization scheme called Power Divergence Subtraction (PDS) and which 
is a generalization of the standard MS scheme based upon dimensional regularization^]. 
Essentially the same method was proposed at the same time by Gegelia[H]. Since then 
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Mehen and Stewart || have made this off-shell scheme (OS) more well-defined and shown 
that it is in fact equivalent to the PDS scheme. 

As a first application the electromagnetic form factor of the deuteron was calculated jj]]. 
Since then the electromagnetic polarizability and the Compton scattering cross section 
for the deuteron have been obtained within the same framework[Q. The inelastic pro- 
cess of radiative neutron capture on protons was first considered by Savage, Scaldeferri 
and WiseB. Based on the same effective field theory also the three- nucleon system and 



neutron-deuteron scattering are now under investigation [10]. To leading order all these 
results depend only on a few known parameters and are similar in structure to the ef- 
fective range approximation in nuclear physics. Even better agreement with data can be 
obtained in higher orders where a priori unknown counterterms appear. When these are 
determined in one process, they can be used predictively in others. 

In the above applications of the effective field theory to low-energy nucleon interac- 
tions there are no complications due to electromagnetic interactions between the nucleons. 
This is obviously not the case for proton-proton scattering which we have set out to re- 
consider within this new framework Jll[|. At sufficiently low energy the Coulomb repulsion 
between protons becomes strong and must be included by non-perturbative methods. It 
is in the same limit that the proton-proton scattering length is determined and a very 
careful analysis is needed to separate the different effects. We have already succeeded in 
calculating the rate for proton-proton fusion into deuterium which is dominated by the 



Coulomb repulsion, using this field-theoretic approach [12]. The same physics is also im- 
portant in other hadronic scattering processes at low energies and in the nuclear bound 
states like 3 He. Similarly, the Coulomb force will dominate in reactions between particles 
of opposite electrical charge at sufficiently low energies. Eventually it forces such systems 
into atomic bound states. 

Experimentally, proton-proton scattering was the first hadronic process studied with 



the help of accelerators and one obtained early on very accurate data] 13]. Landau and 
Smorodinski were the first to construct a formalism in which the Coulomb interaction 
could be separated from the strong interactions ||14|| . This was completed by Bethe in terms 
of a generalized effective range expansion of the low-energy phaseshift for proton-proton 
scattering] 15]. The more model-dependent connection between the measured scattering 



length and the purely hadronic interaction was obtained in the phenomenological analysis 
of Jackson and Blatt []T6|| . Even today this is a main reference to the electromagnetic 
effects involved, something which reflects the absence of a more modern and fundamental 
description of this process at low energies. The effective theory of Kaplan, Savage and 
Wise represents a new and important step in this direction. 

In the next section we present the effective theory and give a short review of the 
applications to proton-neutron and neutron-neutron scattering where Coulomb effects are 
absent. For the proton-proton system these are included in Section 3. We show that 
a perturbative calculation of these effects breaks down at low energies. The scattering 
amplitude depends on the non-relativistic Coulomb propagator modified by the strong 
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interaction and can be calculated non-perturbatively. After PDS regularization we then 
obtain a field-theoretic derivation of the Jackson-Blatt result for the scattering length 
to leading order in the theory. A similar analytical result is also found using instead a 
momentum cutoff as a regulator. Hadronic scattering channels where the Coulomb force 
is attractive is also considered within the same theory. This result is of importance for the 
calculation of energy level shifts in pionic atoms and pionium[|l7]]. Our leading order results 
have been confirmed by Holstein using both standard quantum mechanics and effective 
field theory with an ordinary momentum cutoff as an ultraviolet regulator combined with 
numerical integration^]. He has also considered the more phenomenologically important 
case of several coupled channels. 

In Section 4 we consider effective range corrections by including the dimension-eight 
derivative interaction to first order in perturbation theory. We are then faced with new and 
more divergent integrals involving the Coulomb propagator. By considering the Fourier 
transforms of the Coulomb wavefunctions, we manage to regularize these within the same 
PDS scheme as previously used. Details of this calculation are given in an appendix. The 
resulting hadronic scattering length becomes in better agreement with what one obtains in 
potential models. It also follows that the effective range itself is not affected by Coulomb 
interactions to this order in the effective field theory. 



2 Effective theory for non-relativistic nucleons 

For nucleon momenta smaller than the pion mass, we can integrate out all other fields 
including the pion field. The effective Lagrangian will only involve the nucleon field N T = 
(p, n) and derivatives thereof [fill . It must obey the symmetries we see in strong interactions 
at low energies, i.e. parity, time-reversal and Galilean invariance. Photons are coupled 
to satisfy local gauge invariance and we will assume here that isospin symmetry is not 
broken. The Lagrangian can then be written as a series of local operators with increasing 
dimensions pi] Q JlS|] . In the limit where the energy goes to zero, the interactions of lowest 
dimension dominate. For this case the relevant Lagrangian is thus 

c = Art ^ + Ypj N _ Co (N T nN) • (A^nAO* (l) 

where M is the nucleon mass. The projection operators IL enforces the correct spin and 
isospin quantum numbers in the channels under investigation. More specifically, for spin- 
singlet interactions LTj = o^^Tj/v^ while for spin-triplet interactions IL; = &2&iT2/ 
This theory is now valid below an upper energy A which is set by the pion mass. It is also 
the physical cutoff when the theory is regularized that way. 

The above contact interaction has dimension D = 6 and is thus non-renormalizable. 
It corresponds to a singular delta-function potential and corresponds to the four-nucleon 
vertex in Fig.l. In order to estimate the importance of higher order diagrams like the 
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Figure 1: The four-nucleon vertex corresponds to a delta- function potential. 

bubble correction to the scattering amplitude in Fig. 2, one needs counting rules. If a 
characteristic energy Q flows through the diagram, the propag ator l/(w - Q 2 /2M + ie) 
scales as M/Q 2 since the energy is typically ui ~ Q 2 . For the same reason the phase space 
factor in the loop integration J duid?q/ '(27r) 4 picks up a characteristic factor Q 5 /4irM. The 
estimated magnitude for the bubble diagram in Fig. 2 is thus C^MQ/Air. This will be a 
perturbative correction when CqMQ/Ait < 1 and only a few diagrams will suffice. No real 
or virtual bound states will form in this case and the scattering length have the natural 
size a ~ 1/A. 




Figure 2: One loop correction to the scattering amplitude. 

More interesting is the situation when CqMQ/4tt > 1 and the physics become non- 
perturbative. All diagrams will then contribute to the same order in Q provided the 
coupling constant Co runs in such a way that Co oc 1/Q. Kaplan, Savage and Wise 
showed that a renormalization parameter \/a < fx < A must then be introduced in their 
PDS regularization scheme^] . The scattering length a can now be unnaturally large, i.e. 
a ^> 1/A as it is in the nucleon-nucleon system. The full scattering amplitude T due to 
the strong contact interaction is then the sum over all the chains of bubbles in Fig. 3. It 




Figure 3: The higher order corrections to the scattering amplitude, 
forms a geometric series with the sum 

T{p) = T^W) 
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where 



is the bubble integral and E = p 2 /M is the total center-of-mass energy. The integral is 
linearly divergent but is finite using dimensional regularization in a space of dimension 
d < 2. It has a simple pole when d = 2. In the PDS scheme one subtracts this pole which 
requires the introduction of a dimensionfull parameter fx. Analytically continuing back to 
d = 3 one finds 

Jo(p,M) = -(^)(M + *P) (4) 

The scattering amplitude (§) has then the same structure has the S-wave partial wave 
amplitude 

T{p) = - — (5 

One then recovers the effective range expansion for the phase shift 5 

pcot 5 = -- + ^r p 2 + . . . (6) 
a 2 

in the zero momentum limit when the coupling constant takes the renormalized value 

When only the lowest order coupling parametrized by Cq is included in the effective 
Lagrangian (|l]) we see that the effective range vanishes tq = 0. The small inverse scattering 
length has thus a fine-tuned value given by the difference between two large quantities. 
For example, in proton-neutron scattering we have a pn = — 23.7fm in the spin-singlet 
channel p(I]. Choosing the value fi = m w for the regularization parameter, we obtain Cq = 
—3.54 fm . Physical results should be independent of the exact value of the renormalization 
mass fi as long as 1/a < fx < m n parameter, but it strongly affects the values of the 
coupling constants whose dependence on fi is determined by the renormalization group. 

For neutron-neutron collisions the scattering length is a nn = — 18.4fm, i.e. in magni- 
tude 25% smaller than for the neutron-proton system. A small part of this is the result of 
the proton-neutron mass difference which is a purely kinematic effect. The difference in 
the corresponding coupling constants is however much smaller. Taking again fi = we 
see from (^) that a small relative difference AC® /Co is magnified into a larger difference 
Aa/a ~ afiACo/Co in the scattering lengths. This is the well-known 'amplification effect' 
in more standard nuclear physics where the delta-function potential represented by Co is 
replaced by a potential well of finite extension corresponding to I//i[^J|. The above large 
difference between a pn and a pn thus corresponds to a much smaller difference of less than 
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2% in the coupling constants. This is a dynamical effect due to quark mass differences 
and breaking of isospin invariance by electromagnetic interactions at shorter scales [poj]. 
Recently, Epelbaum and Meissner has shown that these symmetry breaking effects in the 
nucleon scattering lengths are mostly due to the pion mass difference occurring in loops 
in higher orders of chiral perturbation theory ]22]]. 



3 Coulomb corrections to low-energy elastic scattering 

At very low energies it is the repulsive Coulomb force which dominates in proton-proton 
scattering. Increasing the energy, it will still dominate in the forward and backward 
directions while for other scattering angles it is overcome by the strong interaction of very 
short range |L3|]. In our case it is described by the singular potential Cq 5(r). The effects of 
transverse photons are negligible since they couple proportionally to the proton velocity. 



Figure 4: The lowest order Coulomb correction on external legs. 

The lowest order Coulomb correction to proton-proton scattering amplitude in Fig. 
1 is given by the Feynman diagram in Fig. 4. Using the counting rules in the previous 
section, it is expected to give a correction 5T to the result Co from Fig.l of the order 

(M\ 2 e 2 Q 5 aM 

where the characteristic momentum Q is here the proton momentum p in the CM frame. 
This should be compared with a direct calculation of the Feynman diagram which gives 

r d 3 k e 2 M 
J (2tt) 3 k 2 + A 2 p 2 — (p — k) 2 + ie (9) 

where the Coulomb photon has the mass A — > which acts as an infrared regulator. The 
integral is then finite. Combining the two denominators with the Feynman trick and using 

d d q 1 _ A-2-P T(p-d/2) 



we find the result 



(2vr) d (q 2 + A)p (4vr) d / 2 T(p) 



ie 2 f 1 1 
5T = -C M — / dx - (11) 
8vr Jo y'x 2 p 2 — (1 — x)A 2 

= -Cb^f£ + iln^i+0(A) (12) 



2p \2 A 
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It is in agreement with the estimate in (|8|). Since the term which depends on the photon 
mass is imaginary, it will not contribute to the scattering cross section since \Cq + 6T\ 2 = 
Co (Co + 2Re ST) to this order. The cross section is thus infrared finite and proportional to 
1 — nr] where rj = aM/2p. Including one more Coulomb photon exchange in the diagram 
Fig. 4, we find from the same power counting rules that it will contribute a term of the 
order Co(aM / Q) 2 ~ Cor/ 2 . For momenta p < aM/2 we will have rj > 1 and perturbation 
theory is seen to break down. The Coulomb repulsion is then strong and must be included 
in a non-perturbative way. 



Figure 5: The first Coulomb correction on the internal bubble. 

In the same way as the external particles are strongly influenced by the repulsive 
Coulomb potential at very low energies, also the interaction with the strong potential 
Co S(r) is much modified. This can already be seen in the first Coulomb correction to the 
bubble diagram in Fig. 2 given by the two-loop diagram in Fig. 5. It gives a correction 
5Iq to the integral (|j) whose size again can be estimated from the counting rules. Since it 
contains two loops and a Coulomb propagator it should be 

(MQ\ 2 e 2 aM 2 

5I ° * UrJ q 2 = ~4^~ (13) 

More accurately, it is given by the double integral 

f d 3 k r d 3 q M e 2 M 

~ J (2tt) 3 J (2vr) 3 p 2 - q 2 + it k 2 + A 2 p 2 - (k + q) 2 + it ( ' 

It is seen to be infrared finite so that we can take the photon mass A = 0. Since it 
is logarithmic divergent in the ultraviolet, it will have a pole in e = d — 3 when it is 
evaluated using dimensional regularization. This is done in Appendix A where we find 

Here [i is the renormalization mass and Ce = 0.5772. . . is Euler's constant. The diver- 
gent term will be removed by counterterms which modify the coupling constant Cq. The 
prefactor which sets the magnitude of this correction, is seen to be in agreement with the 
above counting argument. In the low-energy limit where the proton momentum p — > we 
see again that this two-loop correction becomes large. In order to have a finite scattering 
length, the logarithmic divergent term must be cancelled by a corresponding logarithmic 
term coming from the exchange of two or more photons in the bubble. These higher order 
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Feynman diagrams are both infrared and ultraviolet finite. Each additional photon ex- 
change is seen from the counting rules to bring in a factor aMj Q to the bubble correction. 
When the external momentum goes to zero, perturbation theory breaks down again and 
the Coulomb-corrected bubble must be calculated in a non-perturbative way. 



3.1 The Coulomb propagator and wavefunctions 

In relative coordinates the two-particle nucleon system is described by the free Hamiltonian 
Hq = p 2 /M where the reduced mass is M/2. The propagator or retarded Green's function 
for this system is 

G { +\E) = 1 (16) 

° E-Ho + ie 

where E = p 2 /M is the total energy. When a complete set of plane wave eigenstates |q) 
is inserted in the numerator, it becomes 

The propagator from an position r to a final position r' in coordinate space is therefore 
(r' |G?o| r) = Gq(E; r', r). It follows that the bubble diagram in Fig. 2 which is numerically 
given by (f|), is just the propagator from zero separation to zero separation, Io(p) = 
G (£;0,0). 

Including the repulsive Coulomb potential Vc = e 2 /4-7rr which acts between the pro- 
tons, the retarded and advanced Green's functions are now 

G^\E) = ^-U (18) 

° E-H Q -V c ±ie 

depending on the boundary conditions which are imposed at infinity. This is equivalent 
to the integral equation Gq = G^ + GqVcGq^. After iteration, is gives the Coulomb 
propagator as an infinite sum of Feynman diagrams with zero, one, two and so on ex- 
changed photons as shown in Fig. 6. 

The Schrodinger equation (H — E) \ tp) = where H = Hq + Vc is the full Hamiltonian, 
has the corresponding incoming and outgoing solutions |*0p )■ They can formally be 
expressed in terms of the free solutions as 

\^) = [l + G ( £ ) V c ]\v) (19) 

This is most easily seen when one uses the equivalent expression \ipp ) = G^ Gq 1 \p) . 
These solutions now have the same normalization as the above plane waves so that 
(tpq \ipp >) = (27r) 3 <5(q — p). Explicit solutions in coordinate space is found from solv- 
ing the Schrodinger equation and can be expressed in terms of confluent hypergeometric 
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Figure 6: Coulomb propagator as an infinite sum. 



or Kummer function M(a, b; x)p5|]. For the repulsive Coulomb potential Vc = ct/r the 
in-state solution with outgoing spherical waves in the future is 



( r ) = e -?*vr(i + ir))M(-ir], 1; ipr -ip-r) e ip ' ! 



(20) 



The corresponding out-state has incoming spherical waves in the distant past and is given 
by the wavefunction 



^H(r) = e -i 7r? ?r(l - irj)M(ir), 1; -ipr -ip-r) e ipr 



(21) 



where rj = aM/2p is the parameter which also appeared in the earlier perturbative calcu- 
lation. The probability to find the two protons at zero separation is thus 



Ci 



l^±)(o)| 2 = e -^r(i + i V ) r(i - irf) 



2tT7] 



e 2-Kri _ I 



(22) 



which is the well-known Sommerfeld factor [p^j pi]. When the relative velocity between 
the particles goes to zero, it becomes exponentially small. At higher velocities rj < 1 and 
the Coulomb repulsion is perturbative. We then recover to lowest order the result 1 — ttt] 
obtained from the Feynman diagram Fig. 4 in the previous section. 

With these Coulomb eigenstates we can now find a more useful expression for the 



Green's functions (18). Since the scattering states form a complete set in the repulsive 



case we consider here, we can instead write it as in (17). Taking the matrix element in 
coordinate space, we then have for the retarded function 



(r'|G, 



(+)i 
c I 



M 



d 3 q vi + V)V4 +) » 



(2nf 



p2 _ q2 _|_ ^ e 



(23) 



In the next section we will see that this propagator gives the main part of the non- 
perturbative Coulomb corrections of the strong scattering amplitude. 
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3.2 Scattering amplitudes and the modified effective range expansion 



Including the strong interaction via the local potential operator Vg, the complete Hamilto- 
nian becomes H = Hq + Vc + Vs- From the full Green's function G^q = 1/(E — H± ie) in 
the presence of both the potentials one can then formally construct incoming and outgoing 
solutions 



(±)\ 



{l + G^{V s + Vc)]\v>) 



(24) 



as in (19). Using now the operator identity A 1 — B 1 = B 1 {B — A)A 1 which impl 
the relation 



ICS 



G 



(±) 
sc 



G 



(±) 
c 



G { c ] V s G 



(±) 
sc 



(25) 



between the Green's functions, we can write the above formal solutions of the coupled 
problem in terms of Coulomb states alone, 



(26) 



n=l 



The scattering amplitude is given by the S-matrix element which is the overlap between 
an incoming state with momentum p and an outgoing state p'. It takes the standard 
form[E5|l 



5(p',p) 



(<7 } l4 +) ) 

(27r) 3 <5(p' - p) - 2mS(E' - E)T(p',p) 



(27) 



With the scattering states (^) the T-matrix element can be written as the sum of two 
parts, T(p', p) = Tc(p', p) + T S c(p', p) where 

(p'l^l4 +) ) 



7b(p',P 

is the pure Coulomb scattering amplitude and 

Tsc (p',P 



(^~ ] \Vs\*l +) ) 



(28) 
(29) 



is the strong scattering amplitude modified by Coulomb corrections [26]. 



Since the full Coulomb wavefunction (r) is known, one can calculate exactly the 
scattering amplitude (p3). The result has the partial wave expansion [pi] 



Tc (p',p) 



47T 

M 



E ( 2i + 1) 



E=0 



2ip 



Pi (cos I 



(30) 



where 6 is the CM scattering angle and = argr(l + £ + irj) is the Coulomb phaseshift. If 
the full scattering amplitude T(p', p) is defined to have the phaseshift ai + 5i, the modified 
strong amplitude (^) will then have the partial wave expansion 

Tsc(p',p) = - J} Y,( 2i + 1 ) e2tae 



E=0 



1 



2ip 



Pt(cos0) 



(31) 
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It should be stressed that the phaseshift 5i is not the same as one would have in the 
absence of the Coulomb interaction. But it can be directly measured from the experimental 
differential cross sections. 

In our case the strong interaction potential Vs = CoS(r) and will only affect 5-wave 
amplitude which we denote by Tgc(p)- If the corresponding phaseshift is called 5, we see 



from (31) that they are related by 



4tt e 2ia ° 

„ {colS -i) = --— (32) 

Because of the strong effects of the Coulomb interaction in the low-energy limit p — > 0, 
the effective range expansion @ does not apply to this phaseshift. This was analyzed in 
detail by Bethe[[l5|] and the result can be written as the generalized expansion pgj] 

Clp (cot 5-i) + aMH(ri) = - — + \r p 2 + . . . (33) 
' a c 2 

Here C 2 is the Sommerfeld repulsion factor (|^) while ac and ro is respectively the S- 
wave Coulomb-modified scattering length and the effective range for the elastic scattering 
process under consideration. In the case of proton-proton scattering where the Coulomb 
potential is repulsive, the function H{rf) is 

H(ri) = Mr)) + ^~- ln(irj) (34) 

LIT] 

where the ^-function is the logarithmic derivative of the T-function. It represents the 
effects of the Coulomb force on the strong interactions at short distances. Using the 
relation 

1 7T 

Imip(ir]) = - — h — coth7T7/ (35) 

we see that lmH(rj) = C 2 /2rj and the imaginary parts cancel out in (|33|). The left-hand 
side will then be real and instead involve the function h(r]) = He^(irj) — In 77 which is more 



suitable for phenomenological analysis [16] 



With this formalism one can extract the physical scattering length ac from the exper- 



imentally measured cross-sections for proton-proton scattering at low energies [16]. One 
then finds a value ac = — 7.82fm and ro = 2.83 fm for the effective range pl|[. While the 
effective range is essentially the same as measured in pn and nn scattering, the scattering 
length is in magnitude less than one half the values found in these processes. This is due 
to the Coulomb effects contained in ac- They can only be removed in a model for both 
the strong and electromagnetic interactions at short distances. Describing these forces 
in a potential model and solving the corresponding Schrodinger equation, Jackson and 
Blatt could isolate a strong scattering length as which is determined by the measured one 
through their relation [16] 



— = — + aM 

as a c 



In— 0.33 

aMro 



(36) 
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With the above value for the effective range tq one then finds as = — 17.0 fm which is very 
close to the ran scattering length. The term —0.33 in the above formula is found to be 
only weakly dependent on the exact form of the strong potential^]. 

This modified effective range expansion obviously applies also to other processes like 
ir + p or 7r + 7r + elastic scattering at low energies where repulsive Coulomb interactions are 
important. The only modification needed is to replace the mass M in ( |33|) with twice the 
reduced mass m = miW2/( m l + where the scattered particles have different masses. 
Similarly, for elastic scattering in channels like ir~p or 7r~7r + where the Coulomb force is 
attractive, it also takes the same form. But it then involves a slightly different function 

H(rj) = j,(irj) + ^~~ H-iv) (37) 

AIT) 



where now r\ = —aM/2p is negative[28|. We will see in the following that both of these 



functions arise naturally in the present theoretical analysis. 



3.3 Coulomb corrections in the repulsive channel 

With the formalism established we can now calculated the Coulomb-modified amplitude 
(p9|) where in our case the strong interaction is represented by the contact potential Vq 
with (p'|Vo|p) = Co- The outgoing scattering state |*p ) can from (|26|) be represented 
in terms of Coulomb states as 

oo 

l*i +) ) = £(4 +) W4 +) > (38) 

n=0 



As a result, we then have for the scattering amplitude 

oo 

Tsc{p) = E \% (G^Vo) n \^ P +) ) (39) 



n=0 

To first order in the strong coupling this is just 



T ^ {P) = /(0/^ ( ^ )|q ' )(q,| ^ o|q)(q| ^ +)) (40) 

= C ^ p >(0)^ p + H0) = C C 2 v e 2 ^ (41) 

after insertions of two complete set of free states. The S*-wave phase shift 2o"o comes from 

the relative phase between the two Coulomb wavefunctions ( |20| ) and (|2l]). In the next 
order of Cq we similarly find the contribution 

T^( P ) = ClC 2 v e 2 ^J o{p ) (42) 

where 



d?k f d?k> ,nM + ) (F)lw _ M f d 3 q V4 +) (0)V4 +) *(Q) 
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is the amplitude G@ (E;0,0) for the protons to propagate from initially zero separation 
and back to zero separation. It can be represented by the one bubble diagram in Fig. 7 
containing the sum of all possible exchanges of static Coulomb photons. 




Figure 7: One bubble diagram with Coulomb corrections. 

Including higher order terms in the expansion (^) we see that they form a geometric 
series with the sum 



(44) 



This result for the scattering amplitude is now to be used in (32) which will give the 
corresponding phaseshift due to the strong interaction. The Coulomb phaseshift <Jq{p) 
and the Sommerfeld factor C 2 ^ are seen to cancel out. We are thus only left with the 
evaluation of (|43|). The two wavefunctions gives the Sommerfeld factor C 2 ^ and thus we 
have 



Hp) = m 



d 3 q 27rr](q) 



(2tt) 3 e 27 "?(9) - 1 p 2 - q 2 + ie 



(45) 



The integral is ultraviolet divergent and must be regularized. Writing Jo = J$ lv + Jq m we 
can isolate the divergent part in 



JS 



div 



-M 



d 3 q 2ir7](q) 1 



(2vr) 3 e 2 ^(e) _ i q 2 

which is independent of the proton momentum p. The remaining integral 



J, 



fin 



M 



d 3 q 2-K7](q) 1 



(46) 



(47) 



(27r) 3 e 27rr ?(<?) -\q 2 p 2 -q 2 + ie 
is finite. It can be done by introducing x = 2irr](q) as a new integration variable and using 



dx 



(e x - l){x 2 + a 2 ) 



In 



2tt 



7T 

a 



1> 



2tt 



(48) 



In our case we will have a = 2nir](p) and the finite part of the Coulomb bubble is just the 
if-function (|34|) already introduced in the problem, 



Ji 



fin 



aM 2 

47T 



H(n) 



(49) 
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Using the full scattering amplitude (Q) in the Coulomb-modified effective range expan- 
sion (33), we see that H{r]) cancels out. The scattering length is thus contained in the 
ultraviolet divergent part of the Coulomb bubble, 



4vr / 1 



1 



a c M \C 



jdiv 
J 



(50) 



There is no contribution to the effective range ro from the Coulomb potential when we 
include only the Cq interaction in the effective theory. 



We evaluate the divergent integral 
3 — d it then becomes 



jdiv 
J 



-M 



a, 



using dimensional regularization. With e 

d-3 2 ^(q) 



(2TT) d Jo 



dqq a 



e 2irr)(q) _ 



(51) 



where VL^ = 2ir d / 2 /T(d/2) is the surface area of the <i-dimensional unit sphere. Again 
introducing x = 2nr](k) as a new integration variable, we find 



div 



-M 



2n d l 2 



2 J T(d/2)(2ir) d 
uM 2 ( u Y r(e)C(e 



airM) 



d-2 



„e-l 



dx- 



4^ \aM^J r(^) 
Here we have introduced Riemann's zeta-function with 

C(6) = -i[l+6ln2vr]+0(6 2 ) 



(52) 



(53) 



From the factor T(e) we get a pole 1/e when e — * 0. In addition there is a PDS pole when 
d — ► 2 from the zeta-function. Since C(e) = £(1 + 2 — d) = 1/(2 — d) + in this limit, it 
takes the form 



4 iv {d 



fiM 1 



(54) 



4vr d-2 

According to the PDS regularization scheme [Q] this contribution should be subtracted 
from the result (52) where then the limit d — > 3 is taken. In this way we are the left with 



div 



j; 



for the divergent part. 



aM l 
4tt 



7 + + 1 



-Ce 



47T 



(55) 



We can now use this result in the expression (|50| ) for the measured proton-proton 
scattering length and obtain 



1 

ac 



Air 
MC 



+ (j,- aM 



1 . Mv/vr 3_ 
- + ln^- + l--C £ 
e aM 2 



(56) 
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The ultraviolet pole 1/e must cancel against counterterms which describe short-distance 
electromagnetic and other isospin-breaking interactions due to quark mass differences [p^j. 
They will modify the coupling constant Co which then takes the renormalized value Co(p). 
It can be used to define a new scattering length 

1 _ 4tt 

0O1) "MC Q,) +/i {b7) 

It is not physical in the sense that it can be measured directly and will thus in general 
depend on the renormalization point fi. Coulomb effects on length scales > 1/fi have 
been removed from it. The coupling constant Cq(^) should be within a few percent of the 
corresponding coupling constants for pn and nn scattering. From (^) we see that a(/j,) is 
related to the physical scattering length ac by 



+ aM 



ac 



ln^ + 1 

aM 



-C E 



(58) 



The result is non-perturbative both in the strong coupling and in the fine structure con- 
stant a which is seen to enter in the combination a In a. This is a consequence of the 
Coulomb force becoming strong at very low energies. Depending on the value of the 
renormalization point /x we see that the Coulomb correction can actually become of the 
same magnitude as the strong interaction. Since the scattering length ac also is negative, 
it can have a very big effect on the size of the hadronic scattering length a(/i). 

Our result for the scattering length is independent of the PDS regularization scheme 
we have used in the above derivation. Instead we can use a simple momentum cutoff A to 
make the divergent integral (|46| ) finite. It then becomes 

it Jo e 27n ?(9) - 1 

Changing integration variable to x = 2irrj, it simplifies to 

aM 2 f°° dx 



jdiv 
J 



(59) 



Ji 



div 



2tt 
MA 
'2^2 



vraM/A x(e a 



1) 



(60) 



+ 



aM 2 

47T 



In- 



2A 

aM 



Ce 



Again it is natural to define a strong and cutoff-dependent scattering length a(A) in terms 
of a coupling constant Co (A) which absorbs the term linear in the cutoff, 



4tt 



+ 



2A 



a(A) MC (A) tt 
In this regularization scheme it is now related to the physical scattering length by 



1 



— + aM 

ac 



hi- 



2A 

aM 



C E 



(61) 



(62) 



a(A) 

which should be compared with (pq). This analytical result is in agreement with what 
Holstein obtained by a numerical integration [lg]. Since in this effective theory the pions 
are integrated out, the magnitude of the cutoff A is set by the pion mass m n . 
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3.4 Elastic scattering in the attractive channel 



We will now consider elastic scattering of two non-relativistic particles with opposite elec- 
tric charge. The full scattering amplitude will again be given by the infinite sum (|39| ) 
where now Gq is the Coulomb propagator in the attractive channel. It will involve 
bound states in addition to the scattering states considered previously in the repulsive 



case. Summing the infinite set of bubble diagrams we find as in (44) for the corresponding 
scattering amplitude 

where now the Coulomb parameter r\ = —aM/2p is negative and Cq is the strong coupling 
constant in this channel. When there are open, strong annihilation channels as in pp, it 
will in general be complex. The full, Coulomb-dressed bubble can be written as the sum 
Jo = Jq + J§ where the first term 

*)=£|^ (64 , 

comes from the bound states. Introducing here the scattering energy E = p 2 /M = 
a 2 M/4r] 2 and the bound state energies E n i = —a 2 M/4n 2 together with the probabil- 
ity ^^(O)! 2 = (aM) 3 1 (87m 3 )5io to find the particles at the origin of the bound state, we 
obtain 

- h , . aM 2 ^ 2r? 2 aM 2 r , , x „ , 

A (p) = -7— E 7 27 2 n = "1— WW + ^ + 2C e\ 65 
47r ^— J n(n z + 7] z ) 47r 

where again Ce is Euler's constant. 

The second term Jq in the attractive Coulomb bubble is due to the scattering states 
and has exactly the same form as (j45| ) in the repulsive case. Since r] is now negative, we 
rewrite the denominator in the Sommerfeld factor as 

1 (66) 



e 2irri{q) _ ]_ ]_ _ e -2-wr]{q) 

This factor can then be split into three parts, Jq = Jq iv + Jq ITI + Jq £W , where the two first 
have the same form as in the repulsive case. In particular, we find 

J (2vr) 3 e~ 27TT i(q) -\q> y ' 



which is therefore given by the result (pop. The 1/e divergence we again absorb into an 
unknown counterterm. Similarly, 

-f in _ f d 3 q -2irrj(q) 1 p 2 

= J (2vr)3 e -2^( 9 ) _ i ^2 p 2 _ q 2 + ie ^ 
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is given by the finite result (|49|). The only new integral appearing here in the attractive 
channel is 



J™ w = M 



d 3 t 



-2irrj(q) 



(2n) 3 p 2 -q 2 + ie 

It is most easily evaluated by dimensional regularization which gives 

aM 2 ~ 



jnew 
J 



2tt 



in ^ +in ^ +i 4 c - 



(69) 



(70) 



It does not have any PDS poles so this is the full result. Adding together these three 
scattering state contributions, we obtain 



<J(\ 



[iM aM 2 

47T 



Ait 



ln ^ + ln( -^-2^ + ^ ) + 1 + ^ 



(71) 



Combining this with the contribution ( p5| ) from the bound states, we finally have for the 
full, attractive Coulomb bubble 



Hp) 



[iM aM 2 



47T 



4tt 



(72) 



From the scattering amplitude (^3[) we can now obtain the physical scattering length 
ac when combined with the corresponding effective range expansion (^) in the attractive 



where ^(77) is the previously defined function (|3 



channel. In analogy with (57) we can define a strong scattering length a(/i) in this channel 
directly given by the strong coupling constant, 



1 



4n 



It can be obtained from the measured scattering length through the relation 



aM 



ac 



In ^§ + 1 
aM 



-C E 



(73) 



(74) 



Since it is not a physical quantity, it depend in general on the renormalization point jj,. 
This result can be obtained directly from the corresponding expression ( |5§| ) in the repulsive 
channel by letting a — > —a in front of the parenthesis. It can also be obtained using a 
momentum cutoff instead of PDS. The integral ( |70D will then be slightly different and 
depend logarithmically on this cutoff. The net result will be as in (p^), again with the 
opposite sign in front of the parenthesis. 

The PDS and the cutoff regularization schemes give slightly different results for the 
the scattering lengths, but it will have no consequence when they are used in other physical 
contexts. For instance, in the attractive channel like in ir + ir~ or ir~p one has hadronic 
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atoms bound by the Coulomb potential. The hydrogen-like energy level spectrum will then 
be perturbed by the strong interaction. With the hadronic contact interactions we consider 
her, we will obtain a shift of the S-states which is easily calculated and proportional to 
the scattering length in lowest order [ 29 1. Since it is only ac and not a(/x) or a (A) which 
have a physical content and can be measured in scattering experiments, it must be the 
one which determines the level shift. This was first shown in a potential model calculation 
by Trueman[3C] and follows also directly using the present effective field theory. This has 
now also been shown by Holstein who has extended the calculation to the more realistic 
situation of several coupled channels [p^]. 



4 Effective-range corrections 



Although the scattering lengths a(/i) or a(A) are regularization scheme dependent and 
cannot be directly measured, one would expect that using our results fl58| ) and fl62| ) for 
them in the case of proton-proton scattering, one would obtained values close to the 
measured values of the strong scattering lengths a pn and a nn . After all, we have just 
obtained an 0(a) correction. However, from Fig. 8 where we plot a pp {ij) as function of 
the renormalization point fj,, we see this the scattering length depends strongly on this 
variable and actually diverges when \x is of the order of the pion mass. The same is the 
case with a pp (A) when the cutoff A increases past m n . This is in contrast to potential 



calculations represented by the Jackson-Blatt formula (36) which permits an almost unique 
determination of the Coulomb-free scattering length and essentially independent of the 
details of the strong potential. One should be able to recover this result by including 
higher order interactions in the effective theory and thereby make the calculation more 
realistic. 

While the leading order interaction term in (JlJ) has dimension D = 6, the next to 
leading order terms in the effective Lagrangian will have dimension D = 8. At non- 
relativistic energies and with only S'-wave interactions there is only one such term, 

£ 2 = ^C 2 (N T V 2 UN) ■ (N T UN)^ + h.c. (75) 

where the operator V = 1/2(V* — i V). It corresponds to a potential V 2 with the matrix 
element 

(q|i> 2 |k) = ^(q 2 +k 2 ) (76) 

Treating this operator perturbatively in the channels with no additional Coulomb inter- 
actions, it was found by Kaplan, Savage and Wise[§|] that the new coupling constant C 2 
is given directly in terms of the effective range tq of the scattering amplitude, 
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Figure 8: Dependence of the inverse scattering length on the renormalization point. Dashed curve 
gives result with only Co interaction, while the solid curve also includes the Ci interaction. 



In the previous section we saw that the leading order coupling constant Co was modified by 
Coulomb effects. Here it will be shown that C2 is unaffected to the order we are working. 



The initial and final scattering states are again given by the interacting states (26) 
constructed from the pure Coulomb states (20) and (^l]). Since we will only be concerned 
with £ = states, it is convenient to use the partial wave expansion 

00 

^(r) = 5^(21 + l)i l Rf\pr)P t {cos6) (78) 
e=o 

which gives for the S'-wave 

R^\pr) = C v e ±ia °e ipr M(l + irj, 2; -2ipr) (79) 

From the Kummer identity M(a, b; x) = e x M(b — a, b; —x) it then follows that the ingoing 
and outgoing wavefunctions are simply related by complex conjugation and differ only by 
the phase factor e iao . 

Treating the new interaction ( |75| ) in first order of perturbation theory, we find the 
Feynman diagrams shown in Fig. 9 contributing to the scattering amplitude. The first 
diagram Fig. 9a gives the amount 

S4%) = \%\ = f /(2^/|r? ^ ] ^ k ^ k2 + k ' 2 H +) ( k ) ( 8 °) 
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Figure 9: The next to leading order diagrams. 
Introducing the quantities 

MP)^ j /^3^ +) (k) = C r) e^ (81) 

and 

MP)^ /(03 k2 4 +) (k) (82) 

we thus have 

dT^ip) = C 2 Mp)Mp) (83) 

The next two chains of bubble diagrams in Fig. 9b and Fig. 9c form similar geometric series. 
Besides the Coulomb-dressed bubble integral (|43|), these diagrams also involve the related 
integral 

j2(p)= /(0/lr¥ k2(k| ^ +)(i?)|k/) (84) 
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Their sum can then be written as 

6T sc C) (P) = ! ^CoMp) Mp)[Mp)Mp) + MP) MP)] (85) 
Similarly, the diagrams in Fig.9d sums up to give 

6T ^ {P) = Y^li^ MP)MP)MP)J2{P) (86) 
Adding up these three partial results, we then have for the perturbed scattering amplitude 

stsc(.p) = {1 ^jll p)] 2 W*(p) + c ° IMpWp) - Mp)Mp)Y) (87) 

Except for ipo(p) which is finite and Jo(p) which has already been evaluated, both V>2(p) and 
J2{p) are divergent and must be regularized. Both of them involve Coulomb wavefunctions 
and there is no obvious way how to do that. 



4.1 Wavefunction regularization 

Since the coupling constants Co(/i) and C2(//) in leading order are known in the PDS 
regularization scheme, it would be simplest if the divergent quantities in (^7j) also could 
be regularized in the same scheme. Since this is defined in momentum space, we will then 
need the Fourier transforms tp^Ck.) of the Coulomb wavefunctions (|78]). These were first 



derived by Podolsky and Pauling pl|. Following them, we show in Appendix B that the 
result can be written as 

OO POO 

^(k) = 4vrV(2^ + l)P*(cos0) / drr 2 Rf ] {pr)j e (kr) 
1=0 Jo 

where is the polar angle in momentum space. As a check, we can then verify that V'o(p) 
as defined in (^) takes its correct value. A similar calculation in Appendix B then gives 
for ip2 (p) the result 



p 2 - naM - ^{aM) 2 



Mp) (89) 



It is obtained using dimensional regularization and the middle, //-dependent term follows 
from a PDS pole in d = 2 dimensions. 

Using the standard representation ( |2"3| ) for the Coulomb propagator, we can now use 
this result to evaluate the double integral (|84| ) for J2. Since the integral over k' equals the 
value of the wavefunction at the origin, we have 

T(n] - M f d " k f d3 i k 2 v4 +) (k)y4 +) *(o) rqm 

MP) - M JJ^f] (2vr)3 - q2 + ie (90) 
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Now we can replace the integral over k by ip2(l) with the regulated result in (|89|) . This 
does not have to be entirely correct, replacing a part of a doubly divergent integral with 
a finite, regulated expression. The problem lies in that ip2(q) has in general higher order 
terms going to zero when e = 3 — d — > 0. However, these can give finite contributions when 
combined with the other divergence in the integral for J2. So we will proceed under the 
assumption that these potential terms are higher order in a so that they can be neglected 
within the accuracy we are working. Thus we have 



Up) = mJ 



(2tt) 3 p 2 — q 2 + ie 



[ p 2 _ mM - i(aM) 2 ] J (p) - J (91) 



where the last term 



J (2vr) 3 e 2 ^(q) - 1 v ; 

is independent of the external momentum p. As shown in Appendix B, it contains a PDS 
pole and will thus be linearly dependent on the renormalization point /i. 



4.2 Effective range and scattering length 

Since we always will choose the renormalization mass /U such that \i 3> aM, we see from 
(jS9|) that if)2 = (p 2 — /iaM)^. As a consequence, the last term V'o^ — 1P2J0 m (|87|) can 
be neglected compared to the first term. Finally we are thus left with the correction 

ST ^ ) = [i-cluvA r2 -'' aU \ <93) 

to the leading order scattering amplitude Tgc(p) m (H)- The p 2 part within the brack- 
ets codes for information about the effective range in proton-proton scattering while the 
momentum-independent part will give corrections to the scattering length. In order to see 
that, we need the modified effective range formula (p3|) which now will appear in the form 



l~ / 1 6Tsc\ 1 1 



-Cy^ «l + aM H( V ) = + -r oP 2 + . . . (94) 

M \J-SC J-SC J a C 1 



Using now 



Tsc C 2 (n) C 2 



p 1 - vaM (95) 



and the ordinary result (|77|) for the coupling constant C*2(//) with the same scattering 
length a = a(/i) which appears in the physical scattering length (|74"|), we see from com- 
paring the 0(p 2 ) terms of this equation that the effective range rg is not affected by the 
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Coulomb interactions to this order in perturbation theory. This result is also in agreement 
with the measured values found for tq in nn, pn and pp reactions p0| f2l|j , 



However, the momentum independent terms in the correction (|93| ) will modify the 



scattering length found in leading order. Instead of (74) we now find 



(/*) 



pp 



+ aM 



1^ + 1 

aM 



;C E 



(96) 



The new terms proportional to ro reduces the downward trend when the inverse scattering 
length is plotted as function of \i as in Fig. 8. When is around the pion mass, we now 
get a value for a(fi) which is much closer to the measured values in pn and nn elastic 
scattering. In fact, we find a(n = m n ) = — 29.9fm which should be compared with the 
value a = — 23.7fm in the pn channel. The differences are now within the range which can 
be explained by electromagnetic interactions at shorter scales [^]. Corrections from higher 
order operators in the effective theory will have a negligible effect on this result since they 
will be associated with higher powers of a. However, they can be slightly different in other 
regularization schemes. 



5 Conclusions 



We have shown that Coulomb effects in proton-proton scattering and other hadronic sys- 
tems at low energies can be calculated systematically in a non-perturbative way based 
directly upon the full Coulomb propagator within the effective field theory of Kaplan, 
Savage and Wise for nucleons. The approach is straightforward and can also be applied in 
other non-relativistic field theories. The phenomenological important quantities in these 
systems are the scattering lengths and effective ranges. While the Coulomb force strongly 
perturbs the scattering length in proton-proton scattering, its effect in irN and tttt reac- 
tions is small due of the small reduced mass in these systems. This is in agreement with 
the experimental situation. Our results are derived in next to leading order and we find no 
changes in the effective ranges in different hadronic reactions due to the Coulomb force. 
The measured values support this conclusion. 

It has been shown that in higher orders in the expansion of the effective theory one is 
faced with increasingly divergent integrals involving the Coulomb propagator. They have 
here been calculated using a method based upon the Fourier transforms of the wavefunc- 
tions combined with PDS regularization. This approach should be put on a firmer basis 
or replaced by a more direct method, perhaps in coordinate space. Also it is of interest 
to do these integrals in other regularization schemes such as the introduction of a simple 
momentum cutoff as we already have used in the simplest case. 
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7 Appendix A 



The twodoop integral ( p~4| ) which gives the lowest Coulomb correction to the bubble dia- 
gram, is infrared finite so we can safely take A — > 0. It is most convenient to evaluate it 
for Euclidean external momentum p 2 = — 7 2 and write it in d = 3 — e dimensions on the 
form 

f d d k r d d q e 2 M 2 

~ J (2ir) d J (27r) d q 2 (k 2 + 7 2 )[(q + k) 2 + 7 2 ] { ' 

after a shift of integration variables. We combine the last two factors with the Feynman 
trick, integrate over k using ( |l0| ) and get 

SIo _ rp-ff) f' dx fJ^J^L (98) 

(4vr) rf / 2 J J {2n) d q 2 A 2- d 



2 



where A = x(l — x)q + 7 • Using now the more general combination formula 

1 f 1 , u a - l uP- 1 , , 

~\77~a TT\ (99) 



A a BP Jo [ujA + (1 - u)B] a +P 
it follows again from the general integral (|l0| ) that 



A i y 1 . !- d/2 r(3-d) (u;a) d - 3 

(2<ir) d g 2 (q 2 + a) 2 - d l 2 Jo T(2 - d/2) (4vr) d / 2 { ' 

where a = 7 2 /x(l — x). Collecting the different factors, we then have for the twodoop 
diagram 

SI = e 2 M 2 ^fc# 7 2 ( d - 3 ) f'dxix-x 2 ) 1 -^ 2 PduJ- 2 (101) 
(47r) d Jo Jo 

The last integral is simply 2/(1 — e) while the first is 



1 , , 2\--4-- f 1 i 1 + e In Va; — a^ 2 



dx(x - x 2 )"2 + 2 = / da; p= = + C(e 2 ) (102) 

o Jo \J X — X 2 
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Now 



dx 



V. 



x — X' 



and 



1 In y/x — x 2 
dx ■ 

'x — x z 



v. 



7T 



-2vr In 2 



We can then take the limit e — > in ( |101| ) and thus obtain 

SIo 



aM 2 f 1 aJ^ „ N 



(103) 



(104) 



(105) 



where Ce = 0.5772 ... is Euler's constant. Finally, we go to back to the physical situation 
by taking 7 = — ip which gives the result (15). 



8 Appendix B 

We will here derive the regulated results (^9|) and (|9l]) from the Fourier transformed 
Coulomb wavefunctions 

V4 ±} (k) = Jd 3 riP^\r)e- ik - r (106) 

with ipp (r) defined by ([ft]). Attaching a hat to the spherical angles of k, we have 

k • r = Ax [cos # cos # + sin 8 sin 9 cos ((f) — <jj)\ (10?) 

For convenience, we choose eft = 0. The integral over the azimuthal angle eft then gives 
simply a Bessel function Jo(—kr sin 8 sin 9) . For the integral over the polar angle 8 we use 
the result 



J d8sin9P e (cos8)J (-krsin8sind)e- lkrcosecose = i e J—Pe{cos8)J e+ i{-kr) (108) 

from Podolsky and Pauling^]. With Jf (—z) = (— X) 1 Jg(z) and introducing the spherical 
Bessel functions 

h{z) = ^J i+ l{z) (109) 

we thus get the result fl88|). 

As a check, we now calculate ipo(p) as defined in (^l|). Using dimensional regulariza- 
tion, it will follow from 

*w=(f)7#* ,(k) (uo) 
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in the limit e = 3 — d — > 0. The angular integration will pick out the £ = part of the 
wavefunction and give a factor An when we take d — ► 3. Thus we are left with a regulated 
expression for the radial integral, 



2 

7T Jo 



<irr Re(pr) / (ifcA; e jo(kr 



The integral over k gives the result 

dkk2~ e Ji(kr) = 22 t 2 



r(l) 



(in) 



(112) 



It is divergent in the limit e — > 0. The remaining integral over r in (|111| ) now involves the 
hypergeometric function (79) and is obtained from the more general formula[23| 

POO 

/ drr^e-^xF^b-pr) = + l)/i"" _1 2-Fi(a,i/ + l,6;p//i) (113) 

Here it gives 

("OO 

/ drr~ 1+e e ipr M(l + ir/, 2; -2tpr) = T(e)(-ip)~ e 2 F 1 {l + iry, e, 2; 2) (114) 

The first factor is seen to cancel the divergence in the previous integral. Combining the 
different terms with 2-^1 (1 + w?j 0, 2; 2) = 1, we finally have 



MP) = 2^e iCTo vr- 1 / 2 r(3/2) 
which indeed is the correct result (|8l"|). 



(115) 



The regularization of ip2(p) in ( p2|) now follows exactly along the same lines and does 
not involve any new integrals. We then find 



rp2(p) 



0F2 e r(-i + 



2-e 



,Fi(l + if7,-2 + e, 2; 2) 



(2k) 



2-d 



K 



d/2 



T(d/2) 



We have here kept the <i-dimensional integration volume element in the last factor. The 
result is now well-behaved in the limit e — > 0. Inserting 2-^1(1 + — 2, 2; 2) = 1/3 — 2r] 2 /3, 
we get 



■02 0)1 



— C»7 e 



p 2 -i(aM) 2 



(116) 



However, this time we see that the above expression also has a PDS pole when d —* 2. Its 
residue involves 2-^1 (1 + ^> — 1> 2; 2) = — ir/ and it thus becomes 



Mp)\^ = C v e^ 2 



(117) 
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It must be subtracted from the previous result with d = 3. We therefore finally have 




which is the result we make use of in (|87|). In the main text we drop the last term since 
we will always choose /i 3> aM. 



The integral (|92|) is of exactly the same form as in (46) and with dimensional regular- 



ization takes the form 




Changing the integration variable to x = 27rr](k) it becomes a standard integral with the 
value 



Making now use of £(— 2) = 0, we see that it is finite when e — > 0. But again there is a 
PDS pole when d — > 2. After having subtracted this part, we get 



Under the assumption fi 3> aM which we make in the main text, this term again can be 
neglected. 
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